Preprint typeset in JHBP style - HYPER VERSION 



MPP-2012-77 



Non-Equilibrium Field Dynamics of an Honest 
Holographic Superconductor 



Xin Gao" '^, Matthias Kaminski^ Hua-Bi Zeng'' ^ Hai-Qing Zhang'^ 

" Max-Planck-Institut fiir Physik (Werner-Heisenberg-Institut) Fdhringer Ring 6, 
80805 Munchen, Germany 

^ Department of Physics, University of Washington Seattle, WA 98195-1560, USA 
Department of Physics, Nanjing University Nanjing 210093, China 
State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese 
Academy of Sciences, P.O. Box 2735, Beijing 100190, China 
'^School of Mathematics and Physics, Bohai University JinZhou 121000, China 
Email: gaoxin@mpp.mpg.de, mski@uw.edu, zeiighbi@gmail . com,hqzhang@itp. ac . cn 



Abstract: Most holographic models of superconducting systems neglect the effects of 
dynamical boundary gauge fields during the process of spontaneous symmetry-breaking. 
Usually a global symmetry gets broken. This yields a superfluid, which then is gauged 
"weakly" afterwards. In this work we build and probe the dynamics of a holographic 
model in which a local boundary symmetry is spontaneously broken instead. We compute 
two-point functions of dynamical non-Abelian gauge fields in the normal and in the bro- 
ken phase, and find non-trivial gapless modes. Our AdSs gravity dual realizes a p-wave 
superconductor in (1+1) dimensions. The ground state of this model also breaks (1 + 1)- 
dimensional parity spontaneously, while the Hamiltonian is parity-invariant. We discuss 
possible implications of our results for a wider class of holographic liquids. 
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1. Introduction 

The gauge/gravity correspondence [1], also known as AdS/CFT correspondence or "holog- 
raphy", has recently provided the tools to study models of strongly correlated systems 
exhibiting superconducting phases [2-4] . All of these constructions have in common that 

first a superfluid emerges from a process which spontaneously breaks a global symmetry. 
Then the resulting field theory is "weakly gauged" afterwards. In other words there are 
no dynamical gauge fields involved in the symmetry-breaking process, which means that 
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these models neglect effects of order 0{1/ gYM^) due to loop corrections which would in- 
volve photons in a real superconductor (^ym is the gauge coupling). This may be seen as 
a shortcoming when trying to model "honest" superconductors in which dynamical gauge 
fields are present during the symmetry-breaking. In this work we address this issue. 

Our first main goal in this paper is to study spontaneous symmetry-breaking in pres- 
ence of dynamical gauge fields. In this sense we are building an "honest" holographic 
superconductor. S-wave superconducting ground states in presence of dynamical gauge 
fields have been studied in higher dimensions [5-8]. In the present work we extend these 
studies in three ways: i.) we examine the dynamics of such a setup, most prominently the 
two-point functions of dynamical gauge-fields, ii.) we consider a vector order parameter, i.e. 
a p-wave solution, and iii.) we work in (1 + 1) dimensions. Note that the Mermin- Wagner 
theorem is evaded in the large N limit [9,10], allowing spontaneous symmetry-breaking in 
two spacetime dimensions. The normal phase of (1 + l)-dimensional fluids in presence of 
an Abelian field was studied holographically in [11], but using an interpretation distinct 
from ours.^ 

Our second main goal in this paper is to study the two-point functions of our dual 
field theory as mentioned in point i.) above. These are in particular interesting since 
our setup forces us to work in the so-called "alternative quantization". We are exploring 
this alternative formulation and its low-energy effective field theory, which we call "non- 
equilibrium field dynamics" . Gauge fields in AdS^, can in general have Neumann or Dirichlet 
boundary conditions. But in this paper we restrict ourselves to a special case where only 
Neumann conditions are possible. We consider a Yang-Mills action without any other 
terms, in particular we have no Chern-Simons term. In this special case gauge fields in AdS^ 
can only have Neumann boundary conditions, so the dual field theory is always written in 
the alternative quantization. This is discussed in greater detail below towards the end of 
this introduction. In this context we also note that we can only study the superconductor 
(dynamical boundary gauge field), not a superfluid (external non-dynamical boundary 
gauge field) in our setup. Therefore we can not directly compare the two within our 
setup. But we find that in presence of our dynamical gauge field a spontaneous symmetry 
breaking occurs in our model. And formally this symmetry breaking in our holographic 
superconductor is very similar to the symmetry breaking known from higher-dimensional 
holographic superfluids. 

We will focus here on the p-wave solution to a Yang-Mills action in AdSa-space 



where g is the metric of AdSs, and F is an SU{2) field strength. Note that we allow no 
Chern-Simons term here. This corresponds to forcing the chiral anomaly in the dual gauge 
theory to zero [12,13]. An action such as (1.1) arises in a particular string theory setup, see 
section 3.1 for details. It can be thought of as a small-field expansion of the non- Abelian 
Dirac-Born-Infeld action of two coincident probe D3-branes in a background generated by 

^We will return to this work for comparison below. 




(1.1) 
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N D3 branes with N ^ oo [16]. For earlier discussions of this brane setup sec [11,17,18].^ 
The Yang-MiUs-theory with action (1.1) is interesting to us for three main reasons: 

First, the expectation is that we introduce gapless modes into an otherwise simple 
setup. For example, in [18] it was found that the duals to AdSa Maxwell-Chern-Simons 
theory do not show any hydrodynamic behavior, all correlators arc fixed by symmetries. 
Now let us naively think of our bulk Yang-Mills-theory defined by (1.1) as such a Maxwell- 
Chern-Simons theory but with vanishing Chcrn-Simons coupling. When we introduce the 
p-wave condensate, wc introduce gapless modes (emerging from the Goldstone modes) 
into an otherwise empty system. So the hope is to study separately these gapless modes, 
possibly analytically. 

Second, our p-wave solution breaks parity in (H-l) dimensions. That means we have a 
parity-breaking ground state from a parity-invariant action. And the effective field theory 
description of the field theory dual to this setup needs to incorporate parity- violating terms. 
In contrast to the Chern-Simons setups in AdSs [24-26] , here the parity- violating transport 
is not related to an anomaly, but to spontaneous symmetry-breaking. 

Third, Yang-Mills-theory in AdS3 is different in that it allows only one particular 
quantization: the bulk gauge fields A{t, x, z) = A{t, x, z) + a{t, x, z) (background and 
fluctuations) correspond to dynamical gauge fields in the boundary field thcorj'^. This 
boundary theory also contains currents, but they are non-dynamical, as they act as sources. 
In other words, the currents and gauge fields have switched their roles compared to the 
ordinary quantization. This claim is based on gauge-invariance and the presence of leading 
logarithmic terms in the boundary expansion of gravity fields. This allows only Neumann 
boundary conditions for the bulk gauge fields. For a careful analysis see [18,27], (sec also 
comments in [28]). One has to carefully repeat the analysis of [29] but for vector fields 
instead of scalar fields. 

For scalars [29] showed that the two possible quantizations are related by a Legendre 
transformation and change of boundary conditions. Also the scalar operator two-point 
function in one quantization is simply the inverse of the scalar operator two-point function 
in the other quantization.'^ A generalization to vector fields in (2+l)-dimensional conformal 
gauge theories (with AdS4 gravity duals) was carried out in [30]: The two quantizations 
are related by an SL{2, Z) transformation in the bulk. On the gravity side, the two bulk 
theories are dual to each other in the sense that S'L(2,Z) is an extension of the traditional 
electric-magnetic duality symmetry. On the field theory side, there is no obvious symmetry 
relation between the two (respectively corresponding) boundary theories. They are in 
general physically incquivalcnt^ , but related to each other: The application of S'L(2,Z) to 
one boundary theory takes us to the other boundary theory. So we can compute obscrvables 
in one quantization if wc know the other. Now in stark contrast to this, wc arc interested in 
a non-conformal theory with gauge fields in (1 + 1) dimensions where only one quantization 

^Closely related to this arc the constructions [19-22]. For a review of the relevant methods see [23]. 
^Note that changing quantizations amounts to changing the boundary conditions and performing a 
Legendre transformation. 

*It has been shown that in the conformally-symmetric case there exists a group-theoretic equivalence 
between the two different quantizations [31]. 
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exists if we do not allow a Chern-Simons term in our action (1.1). So a possible relation 
of conserved current Green's functions to gauge field Green's functions is less clear here, if 
it exists at all.^ 

In our alternative quantization one has to carefully question some of the concepts which 
are by now well established in the normal quantized versions of gauge/gravity. What we 
do not doubt is that gauge/gravity is a correspondence between bulk fields and operators 
on the boundary. In this generality we are going to use gauge/gravity in this work and 
consequently analyze one-point functions and two-point functions of our operators, by 
which we mean our dynamical boundary gauge fields, denoted by 2t(t, x) for background 
fields, and a{t, x) for fluctuations around the background. So in other words A{t, x, z) + 
a{t,x,z) in the bulk corresponds to z) -|- a{t,x,z) on the boundary. In this context 

we are going to make two further assumptions (which are standard in AdS/CFT): 

1. Gauge-invariant bulk field combinations correspond to physical operators in the 
boundary field theory. 

2. Quasinormal modes of gauge-invariant bulk field combinations coincide with the poles 
of two-point functions of physical operators in the boundary field theory. 

Our setup with dynamical gauge fields and fixed currents also obscures the standard 
interpretation in terms of hydrodynamics. Hydrodynamics is a low-energy effective field 
theory which is always expressed in terms of dynamical conserved currents (organized in 
constitutive equations and conservation equations, see e.g. [32] for a review). In contrast to 
this, we have to consider a theory with dynamical gauge fields in the present work. Wc call 
the low-energy effective field theory of our field theory " non-equilibrium field dynamics" , 
or shorter "field dynamics" (in order to distinguish it from hydrodynamics). When doing 
linear response in our setup, we have only access to two-point-functions of gauge fields 
(gauge field propagators), but not to current-current correlators. We note that this formu- 
lation is very similar to the language used in thermal quantum field theory, sec e.g. [33], 
and when discussing non-equilibrium thermal dynamics. Therefore one may hope that the 
gravity setting with Neumann boundary conditions provides results which compare more 
directly with field theory approaches. We stress again that the main point of our work 
is not a comparison of superfiuids with superconductors. Instead we want to explore the 
alternative setup. In particular, we explore what we call field dynamics, i.e. the low-energy 
limit of our alternative description. This is analogous to hydrodynamics in the ordinary 
description. 

We start in section 2 with a few thoughts on possible field theoretic descriptions of 
systems with dynamical gauge fields in two spacetime dimensions. In section 3 we introduce 
a concrete holographic setup, the Yang-Mills action in AdSs and our symmetry-breaking 
solution with vector hair. Then, in section 4 we present our analysis of fluctuations, 
explaining in detail how standard methods need to be adjusted to the situation in AdSa. 

® According to [31], for any tensor field (including scalaxs and vectors) in a conformal theory in any 
dimension the two quantizations are related by a Weyl reflection, and two-point functions are related to 
each other. But our model at hand breaks conformal symmetry. 
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Finally we close with a discussion of open tasks in section 5. Details on gauge-covariant 
and gauge-invariant fields are collected in appendix A. In appendix B we cross-check some 
of our results when using the ordinary interpretation with gauge fields being the sources. 



2. Dynamical Gauge Fields in (1 -|- 1) Dimensions 

Although we are going to study a simple system on the gravity side, there will be several 
complications on the field theory side. We know of no suitable formulation of superfluids 
or superconductors in (1-1-1) dimensions expressed in a way which would be directly 
comparable to our setup. Hence, in this section we will outline our best estimate for such 
a suitable description. In order to disentangle the relevant topics, we separate and discuss 
them here from the field theory perspective (as opposed to the gravity perspective). We 
will be rather schematic and make qualitative observations here, before we introduce a 
concrete model in section 3. 

2.1 Thermodynamics 

As mentioned above, we are going to study gauge fields and their dynamics in this work. 
Finding a version of hydrodynamics with dynamical gauge fields (instead of the usual 
dynamical currents) turns out to be a bit tricky (as discussed below). However, our equi- 
librium states will have a fairly simple standard interpretation. For later convenience we 
here review the canonical and the grand-canonical ensemble. The grand-canonical ensem- 
ble is defined to be in contact with a particle reservoir. Therefore it has a fixed chemical 
potential /v,, and the grand-canonical potential density is il(T, /i). The canonical ensemble 
has a fixed charge density p and the relevant potential is the free energy density F{T,p). 
A Legendre transformation takes us from one ensemble into the other 



So here we have two distinct setups: First, the grand-canonical ensemble (at fixed chemical 
potential p) which we will associate with the "ordinary" setup later. Second, there is the 
canonical ensemble (at fixed charge density p), which we will associate with the "alterna- 
tive" setup later. In this context these two setups have been used many times before, and 
usually one does not talk about this as "two distinct quantizations". But the two setups are 
indeed related by a Legendre transform on the field theory side (plus appropriate choice of 
boundary conditions fixing either fj, ot p). Also the gravity duals to both setups are related 
by a bulk Legendre transformation (plus a change of boundary conditions), which takes us 
from one quantization of the gauge field to the other. This has been first discussed for a 
scalar in [29]. In the stationary state which captures thermodynamics, it is thus clear how 
to interpret the two distinct setups, the "ordinary" and the "alternative" one. 

2.2 Hydrodynamics vs. Field Dynamics 

Let us consider the linear response of a normal fluid at low momentum and frequency in 
1 + d spacetime dimensions 



= F — jip . 



(2.1) 



(j)=G,f -a 



(2.2) 
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where the angular brackets indicate that (j) is a one-point function of the conserved current 
j = 0*! j")"^) a-nd ci = i^t: ^a)^ is the gauge field, which acts as a source here. Greek indices 
run over the d spatial dimensions. By G-^ we mean the (1 + d) x (1 + d) matrix of two-point 
functions of the current j with itself. We define this to be the "ordinary" setup. This is the 
setup which has been realized in most holographic models (including those of holographic 
superconductors) to date. The low-frequency and low-momementum limit of this setup is 
simply ordinary "hydrodynamics". 

Alternatively, we can define a different system by switching the roles of sources and 
one-point functions, so that we get 



where now the angular brackets indicate that (a) is a one-point function of the gauge field a, 
and j is a conserved current, which now acts as source. By we mean the (l + d) x (l + d) 
matrix of two-point functions of the gauge field a with itself. We will refer to this system 
as the "alternative" setup. This is the setup which our gravitational setup (with Neumann 
boundary conditions on the gravity fields) will correspond to through the gauge/gravity 
correspondence. We refer to the low- frequency, low-momentum limit of this setup as "field 
dynamics" (in order to distinguish it from the hydrodynamic description above). 

It is not obvious, how the ordinary setup is related to the alternative one, i.e. how "hy- 
drodynamics" is related to "field dynamics". It is also not a priori obvious what physics is 
captured in the alternative setup at low frequencies and momenta (i.e. in "field dynamics"). 
So one of our goals will be to study the low-frequency, low-momentum fluctuations of our 
gravity model in the alternative setup. As noted before, in our AdS^ gravity model wc have 
no access to the ordinary setup. Therefore a direct comparison between the (well-known) 
ordinary "hydrodynamics" and the (less conventional) alternative "field dynamics" setup 
is not possible in our case. 

Note also that it is not possible to naively invert ordinary hydrodynamics in order to 
get an alternative version of hydrodynamics in terms of gauge fields. In order to see this, 
let us consider a simple ordinary hydrodynamics setup with a constitutive equation for the 
dynamical current 



where = u'^{d^ai, — d^^a^) is the electric field, fi is the chemical potential which couples 
to the charge density p, and = (1,0) is the 2- velocity (all to first order in derivatives). 
Let us assume small field fluctuations a* and a^, to which the system responds by charge 
fluctuations in p. We keep the temperature T and the 2-velocity fixed. Now our goal 
would be to solve this equation for a'^ (such that we could write o'* = . . . , where . . . are 
terms which do not explicitly depend on a'^). Using the conservation equation V^(j'^) = 0, 
and Fourier-transforming by a'^ a'^ exp(—iujt + ikx), we may rewrite (2.4) as a matrix 
equation 



(a) = G«-j, 



(2.3) 




(2.4) 




(2.5) 



-6- 



with the diffusion constant D = cr/x, where x is the electric susceptibility. For this system 
to be invertible, we would need the matrix W to have an inverse. But as we easily see, this 
is not the case since det W = 0. 

2.3 Superfluid Field Dynamics 

In ordinary superffuid hydrodynamics we have to add the phase of the order parameter 
as a new degree of freedom compared to the normal-fluid hydrodynamics. This is due 
to the Goldstone theorem which tells us that for every spontanously broken continuous 
symmetry a massless bosonic field appears in the spectrum of our theory. Superfluids just 
below the critical temperature can be thought of as a mixture of two phases, a superfluid, 
and a normal-fluid one. Once the symmetry-breaking has occured, superconductors can 
be described as charged superfluids for the purpose of studying low-energy excitations 
and transport phenomena. Therefore in the literature superconductors and superfluids are 
often described in the formalism of superfluid hydrodynamics. This is the ordinary setup 
which we will have no access to in our AdSs model. Hence we are not able to formulate 
our results in this common way. 

In our alternative setup we still expect new degrees of freedom since we are still break- 
ing a continuous symmetry spontaneously. To be more precise, our gravity solution will 
develop vector hair below a critical temperature, i.e. we have a non-trivial background 
gauge field Al.{z). This breaks a residual U{1) gauge symmetry spontaneously. Thus we 
expect to see the Goldstone-modes from this symmetry-breaking in the two-point functions 
of our dynamical gauge fields. But it is not clear what the behavior of these modes is going 
to be. For example one could study their dispersion relations uj{k). In the broken phase 
however, we restrict ourselves to the case of vanishing spatial momentum, and compute 
selected Green's functions. Only for the normal phase we consider non-zero momentum 
and find a linear dispersion, as discussed below. 

2.4 Parity- Violating Field Dynamics 

Our gravity solution leads to spontaneous parity-breaking in the dual (1-1-1) dimensional 
field theory. Therefore we also expect to see effects of parity- violation in our low-frequency, 
low-momentum modes. Note that parity- violating hydrodynamics (in the ordinary formula- 
tion) for superfluids with a scalar order parameter has been studied in 3-1-1 dimensions [34]. 
There, parity-breaking leads to a large number of additional transport coefficients. In prin- 
ciple it should be possible to carry out a similar analysis in the alternative setup in order 
to derive a formulation of parity-violating field dynamics. For example one could use the 
fluid/gravity framework [35] and compute the constitutive relations (presumably for dy- 
namical gauge fields) and corresponding coefficients explicitly. It is not clear what to expect 
from such an analysis. We note that adding a Chern-Simons term to our action (1.1) would 
lead to a chiral anomaly on the field theory side. This would break parity as well and such 
setups have been examined recently in [18,28,36]. However, this would also change the 
boundary behavior of the gauge fields. 
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Figure 1: In (1 + 1) dimensions the anomaly diagram on the left is identical to the self-energy 
diagram for the gauge field on the right. 

2.5 Massive Gauge Bosons and the Schwinger model 

One may speculate that the physics of dynamical gauge bosons should be best captured 
by thermal quantum field theory, at least at weak coupling. In this case the correlators 
between gauge bosons a discussed around (2.3) should be interpreted as gauge boson 
propagators. Due to the finite temperature we naively expect our gauge bosons to have a 
thermal mass. Furthermore, they should acquire a mass through a Higgs mechanism due 
to the spontaneous symmetry-breaking. The latter mass contribution should depend on 
the condensate which breaks the symmetry. So in particular this mass contribution should 
vanish as we approach the phase transition from the broken phase. But we will see below 
that neither the thermal mass nor the mass from symmetry-breaking can appear in our 
setup. 

Already at zero temperature gauge bosons in (1 -|- 1) dimensions are different from 
higher-dimensional cases. Schwinger discussed quantum electrodynamics in two spacetime 
dimensions [37]. Among other things he found that photons in (1 + 1) dimensions can have 
a mass term which does not break gauge invariance. Furthermore the self-energy diagram 
of gauge fields in (1 + 1) dimensions is identical to the diagram for the chiral anomaly in 
(1 + 1) dimensions, see figure 1. The mathematical reason for this is that there are two 
7-matrices in (1 + 1) dimensions, 7*^, 7^, plus the 7^ = 7^7^, and they obey the relation 
7/^7^ = —^t^i^j^^ with the totally antisymmetric symbol in two spacetime dimensions e^'^. 
Thus in the fermion loop in figure 1 on the left we can replace the coupling to an external 
axial current i7^75 by the coupling to an external gauge field —i'ju- This implies a strong 
constraint: If a gauge field in (1 + 1) dimensions has a mass, then the theory has a chiral 
anomaly, and vise versa. 

We will use this result below in order to interpret the (1 + l)-dimensional field theory 
which is dual to our AdSs setup. Note that for a massive photon in (1 + 1) dimensions the 
Ward identity requires a pole in the gauge boson propagator at two-momentum k"^ = [37]. 

3. Holographic Setup and Symmetry-Breaking Solution 

While the previous section was very speculative, in the present section we are going to 
work with a concrete gravity model. On the dual field theory side, this will realize a p- 
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wave superconductor in (1 + 1) dimensions with dynamical electromagnetic fields, i.e. the 
alternative setup described above. 



3.1 Dirac Born Infeld Action for Coincident D3 Probe Branes 



In the context of probe brane constructions [16] -apart from the DBI-action- the effective 
brane world volume action also contains a Wess-Zumino term. This leads to the Chcrn- 
Simons coupling of the gauge field. However, for particular brane intersections such a 
Chern-Simons term can be absent and only the Maxwell term remains. We consider a 
defect brane configuration which realizes this very situation. This is the D3 — D3' defect 
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Table 1: D'S — D'3' brane configuration which provides only a Maxwell but no Chern-Simons 
term. We consider N ^ oo D3 branes, and Nf = 2 D3' probe branes. 

configuration illustrated in Table 1 . The metric generated by the N background D3 branes 
is the standard AdS^ x (black brane) metric 



H-^/\-f{r)de + dx"" + dy"" + dz"") + H^^^ + r^dUl) , (3.1) 



fir) 



where /(r) = 1-^ , H = {L/A)^ and dO?^ = dO^ + cos^ 9d^'^ + sin^ OdSl Transforming 
the radial coordinate to n = rg/r, the metric becomes 

'rQ\^ -f{u)dt^ ^dx"^ + dy^ + dz^ L'^dv? 
1. 



ds'^ 



D3 



(3.2) 



v?f{u) 

We consider Nf coincident probe D3' branes with a trivial embedding, i.e. 9{r) = 0. Then 
the metric induced on the world volume of the D3' branes is 

-f{u)dt^ + dx'^ 



ds\)^i 



+ 



L^du^ 



(3.3) 



v?f{u) 

Prom this metric we can see that the probe D3' branes cover AdS^ as well as an S\ cycle 
inside the five-sphere. The existence of such AdS^ x embeddings for the D3' branes 
has been demonstrated in [17]. Furthermore we switch on a U{Nf) gauge field living on 
our (1 + l)-dimensional defect. In fact our gauge field lives on the world volume of the 
D3' branes. These probe branes do not wrap any cycle with Ramond-Ramond fiux and so 
the Wess-Zumino part of the brane action does not give rise to a Chern-Simons coupling 
for our U{Nf) gauge field. This argument was first made in [18] (see also [11]) and works 
for Abelian as well as non-Abelian gauge fields. Then up to some constant the action for 
the probe £)3'-branes will be the non-Abelian DBI-action [38], which after exploiting our 
symmetries^ amounts to 

Sd3' = -TozNfStr J dtdxdzd^^/ -det{gab + (27ra')i^a6) , (3.4) 



®See section 3 of [21] for details of the analogous computation for probe D7 branes. 
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with F = dA + A f\ A, and the symmetrized trace Str is taken over U{Nf) representation 
matrices (we have renamed the radial coordinate u ^ z for convenience). Here we choose 
the metric and gauge field to be indepenent of the 5'^-coordinate ^. After expanding the 
square root in (3.4) in small field strengths F the leading contribution is extremized by the 
AdSs black hole background. The subleading contribution merely has the form of the Yang- 
Mills action which we already advertised in equation (1.1). It should be noted here that 
the non-Abelian DBI action is only valid up to fourth order in field strengths [14,15], and it 
has been shown to disagree with corresponding string scattering amplitude computations 
beyond this order. However, in the present paper we are only interested in the leading 
(quadratic in field strengths) order. 

The field theory which is holographically dual to this gravity setup is known and has 
originally been studied in [17], see also [11,18]. This dual field theory is U{Nf ) x U(N)M = 
4 SYM theory coupled to a bifundamcntal hypermultiplet along the (1 + l)-dimcnsional 
defect. In particular there is a dynamical gauge field 2t living on the (1 + l)-dimensional 
defect. This dynamical gauge field is dual to the dynamical bulk U(Nf ) gauge field A on 
the stack of Nf probe branes. Therefore the setup presented in this subsection provides an 
example of the dynamical boundary gauge fields which wc want to study in this work. It 
should be noted that our dynamical boundary gauge field 21 (dual to the bulk gauge field 
A) is still a gauge-independent operator under the original "color" group U{N — t- oo). It 
is only dynamical and hence gauge-dependent under the additional "flavor" U{Nf) which 
we have introduced. 

3.2 Yang-Mills- Action in AdSs 

In the previous subsection we have illustrated how the Yang-Mills action (1.1) may arise 
from a certain string construction. In this subsection we focus on this Yang-Mills action for 
a non-Abelian gauge field A in AdS^. The neutral AdS3 black hole background in Poincare 
coordinates is given by 

7-2 J 2 

ds^ = —{~f{z)dt'' + dx^ + -^), (3.5) 

in which f{z) = l — z^,z = r+/r, and r+ is the horizon of the black hole. The black hole 
temperature is given by 

T=r±. (3.6) 

We are free to set L = 1 and r+ = 1. To this background, we add SU{2) gauge 
fields ^^j, with a = 1,2,3. Wc split this bulk gauge field yl^ in background fields A^, and 
fluctuations around this background. Our Ansatz for the background fields is 

A = ^iz)dtT^ + W{z)dxT^ , (3.7) 

where r" arc the generators of the SU (2). For our fluctuations wc choose the gauge = 
but leave them otherwise unrestricted for now. The field strength is F^^^ = d^A^ — d^A^ + 
gvM^^'^'^-y^n-A.I,. For simplicity we set qym = 1 in the remainder of this work. Note that 
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the value of gvM would be important if we were to relax the probe limit, comparing the 
gravitational coupling to the gauge coupling. 

The equations of motion for the p-wave background fields $ and W are 

= d,{zfd,W) + j^^W , (3.8) 

= d,{zd,^) - j^W"^ . (3.9) 

3.3 Alternative Boundary Conditions 

The boundary condition at the horizon is $(2; = 1) = 0, which is to make ^dt well-defined 
at 2; = 1, while W should be regular at the horizon. The horizon expansion is 

$ = -$i(^ - 1) + ^ (2^1 + Wo^^i) {z - If + 0{{z - If) , (3.10) 

W = Wo + 0{iz^lf), (3.11) 

This expansion has the same structure as in the AdS5/CFT4 and AdS4/CFT3 cases. In 
contrast to that, at the boundary, the story is different. The asymptotic behavior of W 
and ^ at the boundary 2; = is 

W = -a^lnz + (2t^), 

$ = -;jtln2; + (2tt), (3.12) 

where Zt ='■ P and Zx are the sources, and (2lt) =: ji, as well as (21^:) are the vacuum 
expectation values^ of the operator dual to the bulk background field A^. Note that we 
would get a different boundary behavior for the gauge fields had we allowed a Chern-Simons 
term in our action (1.1). In that case the asymptotics depend on the Chern-Simons level, 
see for example [18]. 

In (3.12) we now encounter the alternative quantization mentioned above, which cor- 
responds to the choice of Neumann boundary conditions [5-8]. In AdS4 for example, we 
would be free to interpret the leading term of $ as the source (chemical potential), and the 
sublcading term as the expectation value of the dual current operator (charge density), or 
the other way around. However, this is not possible in AdSa [18,27]. We have to choose 
the Neumann boundary condition. One may be concerned that [18, 27] consider Abelian 
gauge fields for the most part. But it has been noted in [27] and earlier in [30] that the 
asymptotics and hence also the boundary conditions are determined by the linear theory. 
Our non-Abelian terms contribute to non-linear (interaction) terms. ^ Therefore our non- 
Abelian gauge fields in AdS^ should only satisfy the Neumann boundary condition, just 
like their Abelian counterparts. 

So compared to the higher dimensional cases, we observe two changes: First, the 
constant term is not the leading contribution anymore. A logarithmic term appears and 

'^Note that we have supressed gauge indices a,h, . . . since according to our Ansatz (3.7) in the background 
the a = 3 index always appears with the spacetime index /it = i and a = 1 appears always with fj, = x. 

®As a concrete example, boundary conditions have been discussed in [27] for 50(8) non-Abelian gauge 
fields of AdSi supergravity. 



- 11 - 



contributes to leading order. Second, as explained in [18,27], the leading (now logarithmic) 
term has to be interpreted as a source and corresponds to a (non-dynamical) boundary 
current. Therefore, the source in this system is the boundary current -3. The subleading 
(constant) term is the expectation value of the operator dual to the field A. We denote this 
operator by 21. According to the analysis in [18,27], this operator 21 has to be interpreted 
as a gauge field living on the AdSs-boundary. So in summary, our bulk gauge field A is 
dual to a boundary gauge field 21. This boundary gauge field 21 is sourced by the current 

Thus our fixed charge density p = ^t sources the dynamical gauge field ^ = (2tt). This 
particular component (2lj) can be interpreted as a chemical potential (which does not act 
as a source, but is determined by the equations of motion). Meanwhile (21^;) is the order 
parameter of the symmetry-breaking, while Zx is the source term which we demand to 
vanish, for the breaking to be spontaneous. Our operator 21^^ is the spatial component of 
a dynamical boundary gauge field, which acquires an expectation value. In the context of 
thermodynamics our alternative setup has a simple interpretation: since we are working 
at fixed charge density, this is the canonical ensemble (provided we add an appropriate 
boundary term to the action (1.1) as discussed below). 



3.4 Symmetry-Breaking Solution and Free Energy 

For simplicity we choose to work in the probe limit, i.e. the gauge field does not backreact 
on the AdSs metric. After solving the equations of motion (3.8) with boundary conditions 
(3.10), the order parameter (21a;) is plotted versus the temperature in figure 2, both in units 
of Tc- Note that in the canonical ensemble Tc is determined by the only other dimensionful 
quantity in the problem, namely /?, see [3,4]. So we actually plot (21a;) rescaled by (plnpc) 
versus the temperature as measured by T/Tc = Pc/p- 
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Figure 2: Condensate {0)/Tc = {^x) / {p^npc) versus temperature T/Tc 
temperature is given by Tc = 0.046p. 
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Near the critical temperature Tc, the condensate (21a;) behaves hke 

(21^) «i 25.5515rcVl-T/Tc. (3.13) 

The critical exponent is 1/2 which is identical to the one in mean field theory. Note that this 
solution also breaJiS parity in (1 + 1) dimensions. In any odd number of spatial dimensions 
parity is defined as inversion of only one of the spatial coordinates. In our setup this would 
be inversion of the one spatial coordinate x. Our order parameter (21^:) is a vector pointing 
along this x-direction (or against it depending on the sign conventions). Therefore parity 
symmetry, the invariance of our system under inversion of the x-coordinate, is broken 
spontaneously. 

In the following we will compute the free energy difference between the normal state and 
the broken state in order to verify that the broken phase is preferred at low temperatures. 
With the Ansatz (3.7), we can rewrite 

= _ l^^a^Fa^. ^ -zf{d,Wf + j^-'W^ + z{d,Wf . (3.14) 

Then the on-shell action 
S^^ = J ^/^C(fx = 2lim J dtdx {-zfWd^W + z$5^$)^^, - 2 J dzdtdx j^'^W'^ , 

= 2 J dtdx {-lip + p2 log(z)) - 2 J dzdtdx ^^^W^ , 

(3.15) 

can be obtained by integrating our action (1.1) by parts, and using the equations of mo- 
tion (3.8) in the form 

- j dz[zf{d,Wf] = -zfWd,WU=, - j dzj^^W^ , (3.16) 

and 

J dz[z{d,<^f] = z<^d,^^=e - J dzj^^W^ . (3.17) 

In order to place the dual field theory in the canonical ensemble we need to add a boundary 
term [18, 27] 

5bdy = dtdx (V=5^M^"').=o = 2^ dtdxi-zfW'W + z^'^),=o , (3.18) 

where the prime denotes a derivative with respect to the radial coordinate z. Note that 
this contribution (3.18) is identical to a part of the on-shell action (3.15). Now we have to 
add a counterterm 

^counter = 4 J dt dx ^/^F.^F''' log(z) = J dt dx \og{z) , (3.19) 

in order to cancel the logarithmic divergences in (3.15) and (3.18). Here h is the metric 
induced on the boundary. Hence we obtain the renormalized on-shell action 

.5ren = ^os + S^dy + ^counter = -4 J dt J dxpil -2 j dx^ j^'^W'^ , (3.20) 



- 13 - 



where we have used that the source has to vanish, i.e. = is required, for sponta- 
neous symmetry-breaking. Finally the free energy density T is given by the Euchdean 
continuation of 5ren as 



-TSr, 



Afip + 2 dz 



f 



(3.21) 



Note that the free energy of the normal state in our conventions is always zero, because 
in that phase the general gauge field solution is A'l = —plog{z) + Cq and the constant Co 
has to vanish for j4f(z) to be well-defined at the horizon. This situation would change if 
we were to consider backreaction; then we would have finite p and /i. 

The free energy difference between the broken state and the normal state is plotted in 
figure 3, from which we can see that the state with a condensation is thermodynamically 
preferred when T < Tc- The condensate vanishes with a square root behavior. 

Thus our first result is that there is a second order phase transition for the Yang-Mills 
theory (1.1) in the AdSs black hole background at the critical temperature T^. Our broken 
phase develops a (parity-breaking) vector condensate, which corresponds to vector hair of 
the AdSs black hole. 



Ar 




Figure 3: The free energy difference between the normal and broken phase. 



The p-wave solution we have presented here is different from the one discussed in [39] . 
Our lowest energy solution for the field W{z) has one node at some z = z*, while the p- 
wave in AdS4 has no nodes. One may argue that more nodes should correspond to a higher 
energy for the solution [39]. According to that reasoning the zero node solution should 
always be lowest in energy. However, in our AdSs case the zero node solution coincides 
with the solution which has W^(z) = 0. Or in other words: there is no symmetry-breaking 
solution with zero nodes in our setup. 
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4. Fluctuations 



In this work we will not be able to address all the field theory features presented in section 2 
at once. Instead we are going to pick selected modes in order to investigate the dynamics of 
the full system in the normal phase and in the phase with spontaneously broken symmetry, 
i.e. the broken phase. Our goal is to show that there are non-trivial dynamical low-energy 
modes near the phase transition in the broken phase of this system. 

4.1 Fluctuation Equations of Motion 

We start by writing down the full set of equations of motion after the gauge choice = 
is made 

n _ „i" , fecjg^ - {ikal + Wal)^ 

= «x + (^^ + yj + f2 + p ' 

= af + -af 7 af ^ -. — , 

^ / / 

2// ,1 f\ 2' ^2-^^-2 2 kua^ + i{ua^W + {kaj + 2ujal)^) 



= + (- + j)^x + J2 «x + p 

3// 1 3/ k^ + W\^ kujal + W{-ii2ka^t+ujal) + 2ai^) 
U — df n~ 



= + (- + L)al + —4 + ^2 ^ ■ (4-1) 

Variation of the action with respect to the fields yields the following constraints 

= -iua]' - $af - ikfa]! + a]^' , 

= ^a]' - itoaf - f {ikal! - Wa^J + alw'^ - , 

= -iual' - f {WaX + ik4 - alW) . (4.2) 

Expanding the fluctuations (4.1) around the AdS-boundary, we get 

al = a«(«°"'-^'') log(l/z) + a?(^"^) + . . . , (4.3) 
al = a:^^^-^^^^ log(l/z) + aS(^^^) + . . . , (4.4) 

where a = 1,2,3 is the gauge index, (source) denotes the source, and (vev) denotes the 
vacuum expectation value. We follow the same logic which was applied for our background 
fields A below (3.12). So the coefficient of the logarithmic term in (4.3) is interpreted as 
the boundary theory current = gbi^^^^^^^ which acts as a source. The constant term is 

interpreted as the vacuum expectation value of a boundary theory gauge field (o^) = a^^^^^^ . 
Near the horizon these fluctuations behave as 

= (1 - ( 0{1 - z)) , (4.5) 

< = (l-z)-^t(<(«)+0(l_^)) . (4.6) 
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All higher coefficients depend recursively on the three initial horizon values a^^^\ cLx^^^ 
and a^^^^ of the spatial components. 

4.2 Analytical Results 

Gauge-invariant field combinations Above we have constructed a ground state which 
breaks the SU{2) gauge symmetry explicitly to ^7(1)3, which is broken spontaneously by 
the same ground state. Therefore our ground state does not have any gauge symmetry left. 
However, the Lagrangian of our theory still enjoys the full SU (2) gauge symmetry. Thus, 
while the S'J7(2)-symmetry is completely broken for our background fields which describe 
the ground state, the fluctuations still transform under a full SU{2) gauge symmetry in 
the following way 

Sa"^ = 5A« + e^^'M^A'^ , (4.7) 

with the transformation parameters A"(t, x) where a = 1, 2, 3. These A" do not depend on 
the radial coordinate z because we have chosen to work in the axial gauge Az = earlier. 

First, we note that there are three pure gauge solutions for the fields (aj^, a^, a|, a^, af , a^), 
which in this notation read 



(I) 

(11) 
(III) 



( -ioj, ik, 0, 0, 0) , 

( 0, -ioj, ik, 0, W), (4.8) 

( 0, 0, 0, -W, -iu, ik) . 



We have labeled these solutions (/), (//), (III) for later convenience. 
We find the following gauge-invariant bulk field combination: 

4 = 4 + -4 + w . (4.9) 

At vanishing momentum this expression reduces to the gauge-invariant combination con- 
sidered by Gubser &; Pufu in [39] . 

Of course we can build more gauge-invariant combinations such as 

4 = 4 + ^(^^ - ^')«? - ^4 • 

These are discussed in appendix A. 

In the most general case below Tc, all six fluctuations are coupled through the conden- 
sates, momenta and frequency terms in the equations of motion (4.1). Therefore all the 
two-point functions will have the same poles. For a more detailed discussion see [40,41]. 
We will be interested in finding these common poles in order to identify low-energy modes. 
These are the poles at low frequency and small momentum which govern the long-time 
behavior of our system. By one of our inital assumptions we identify the poles in our two- 
point functions with the quasinormal modes of the gauge-invariant bulk field combinations 
defined above. Our quasinormal modes can in turn be found by a method which we 
outline below. In the same way in which all correlators contain the same poles, the bulk 
fields will all have identical quasinormal modes in the fully coupled system of six equations 
of motion (4.1). 
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Gauge- covariant field combinations In the normal phase the symmetry is only broken 
partly by the ground state. There is a U (1)3 remaining intact. Therefore the gauge 
transformation for the fluctuations in this case amounts to 



(4.10) 



where is the gauge transformation parameter of the remaining U (1)3 in the background. 

It is not possible to form gauge-invariant combinations of the fluctuations as we 
did before in the broken phase. However, we can form gaugc-covariant combinations, i.e. 
combinations which will have a well-defined homogeneous transformation under the gauge- 
symmetries. 

We consider the following gauge-covariant field combinations 



63 



($(z) - w)e+ - kef 
{^{z) + uj)e- + ke^ 
ka^ + ijja^ , 



(4.11) 
(4.12) 
(4.13) 



where = aj ziz iaf , = it ia^ . Note that 63 is gauge- invariant and thus corresponds 
to a physically meaningful operator on the boundary by our initial assumptions. Details 
are discussed in appendix A, see also [42]. 

Determinant Method & Quasinormal Modes Quasinormal modes are those bulk 
field solutions at a particular = LOn,QNM which are infalling at the horizon and satisfy a 
Dirichlet condition at the AdS-boundary; we choose them to vanish at the AdS-boundary. 
In a coupled system it is not a trivial task to identify the relevant field combinations which 
have to vanish for a quasinormal mode to be present. A systematic method has been 
developed in [40,41]. However, we will see that this method needs to be extended in our 
case at hand, because there are logarithmic sources. 

Let us solve the equations of motion numerically, for example with the horizon values 
l^^^ = -1, al^^^ = 1, and a^^"^ = 1, see (4.5). We label this solution {IV). Then we 

" = -1, 



solve the equations of motion again using different horizon values a^f^^ = 1, a^^^^ 



and a^^^^ = 1- This solution is linearly independent from the first one and we label it (V). 
There is a third linearly independent solution which we find using the values a^.^^^ = 1, 
a^^^^ = 1, and a^*'^'* = —1, and we label it {VI). It was shown in [40,41], that these three 
solutions can be combined in a matrix together with the three pure gauge solutions given 
in (4.8). The determinant of this matrix, evaluated at the boundary z = was shown to 
yield locations oj^qliMi^) °f quasinormal modes of the coupled system. Quasinormal modes 
can be found by considering the following determinant 



V = det 



—iuj 


ik 


















—ioj 


ik 





w 











-W 


—ioj 


ik 



\ 



h 

l{VI) 



x 

l{V) 

X 
X 



h 

■2{VI) 



2(IV) 

X 

2(V) 

2{VI) 
x 



H 



X 

3(VI) 

X 



(4.14) 
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In [40, 41] this kind of determinant was required to vanish in the following way: 



= lim P . 



(4.15) 



This procedure essentially sets the relevant sources to zero if we work with bulk fields which 
have a constant as the source term, and no more divergent terms. However, in our AdSa 
case we have logarithms. Therefore the recipe has to be changed in order to extract the 
sources which now hide in the logarithmic terms. One sensible prescription appears to be 



Whenever we set a fluctuation consistently to zero, we should eliminate the corresponding 
column and one row of solutions from this determinant. For example this prescription yields 
exactly the gauge- invariant combination defined in (4.9), when we set = = at 
W,^0, but k = 0. 

Let us work through this example: we have two pure gauge solutions a = {—ico, 0) 
and a = {—^,—ico,W), and one non-trivial solution, which depends on one parameter, 
namely the source of at the horizon. So this gives a solution matrix 



Now we compute the determinant V = det H = a^ — VF(<&a| +iwa|) / (w^ — $^), and require 
it to vanish in general. This yields exactly the gauge-invariant combination which was 
found in the analogous AdS4-system [39] . But there is a crucial difference now. In AdS4 we 
simply had to take the prescription (4.15) because the sources were the asymptotic values 
of the fields, and was just a linear combination of these. Luckily the latter is still true 
in AdSs, but now we have to extract the logarithm coefficients. This leads to the adjusted 
prescription (4.16). The reason why this works here is that for k = we can reduce (4.1) 
to three second order equations of motion for aj, af, and a^. These three equations are 
coupled. But gauge symmetry restricts them, so that they can be reduced to one single 
equation. In other words: When we work out all of the 9 equations of motion (not choosing 
a gauge yet), then we see that there are 2 constraint equations from (4.2) which relate aj 
to af and a^. Then the counting is: 3 equations of motion with each 2 boundary conditions 
give 6 boundary conditions. Choosing the infalling solution for each of the 3 fields there 
are 3 boundary conditions left to be fixed. Now we can use our 2 contraints, which leaves 
merely one solitary boundary condition. This means that there is only one relevant field 
and our one boundary condition determines the normalization of this field. 
As a second example, we consider the normal phase with W = 0. Here the adjusted 
method (4.16) yields accordingly 63 = ua^ + ka'l at a^ ,j, = = Oj ^. and $ / 0, and correctly 
requires the leading logarithmic term in 63 to vanish. However, for the genuinely coupled 
cases the prescription (4.16) may have to be modified.^ In particular the gauge-covariant 

^Equivalently, at least in the normal phase, one can carry out an analysis formally similar to [11]. 



= lim zdzV . 



(4.16) 



H 



I -iuj $ \ 
— $ —iuj W 



(4.17) 
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gauge fields will not necessarily be linear combinations of the individual fields anymore. 
We postpone their consideration to future work. 

Numerically we can not take the limit z — > exactly. Therefore we have to choose the 
same cut-off A = 10~^ which we had chosen when computing our background data $ and 
W. See [40, 41] for a more detailed discussion. 

4.3 Numerical Results 

Here we will simplify our lifcs and consider only particular gauge-invariant field combina- 
tions which decouple from all others. Such field combinations are completely determined 
by requiring the infalling boundary condition at the horizon, and by giving a single second 
boundary condition fixing the normalization of the solution. For a discussion of how to 
find these decoupled solutions see section 4.2 and the appendix A. Note that the frequency 
uj in the following plots should be understood as the dimensionless quantity ui ^ uj/ (27rT) , 
and similarly k k/{2TrT). 

Normal phase at vanishing spatial momentum k = We solve the equations of 

motion (4.1) in the normal phase at finite chemical potential induced by $ 7^ at w 7^ 0, 
but at W = and k = 0. In this case the last equation in (4.1) can be consistently 
decoupled from the others. Therefore we consider that single equation for the fluctuation 
a^, setting all other fluctuations to zero 

„ 1 — 0/ o / . _ „s 

Note that this equation does not depend on the background field $ in other words the dual 
operator is not charged under f7(l)3. We consider the Green's function of the dynamical 
gauge field combination in terms of the gauge- invariant electric field ei^ = wa^. The 
straightforward recipe for correlation functions (up to contact terms) can be applied here 
as shown in [18], to yield 

9 (vev) r,o?^ — \ov{z)zdzot 

Gfs^s = oj^G^s^s = oj^j^ = oj^^ ^ 3 " ^ . 4.19 

'^x'^x (source) zOxOg 

Note that we are consequently ignoring numerical factors for the sake of a simple presen- 
tation here. Solving the equation of motion (4.18), we see that the imaginary part of the 
Green's function here asymptotes to a constant at large frequencies and has a pole at a; = 
independent of charge density. 

Normal phase at non-vanishing spatial momentum k ^ We solve the equations 
of motion (4.1) in the normal phase with = at finite momentum k ^ and at finite 
chemical potential induced by $ 7^ 0, while oj ^ 0. Now the last two equations in (4.1) 
decouple from the others. Therefore we consider only the two coupled equations for the two 
fluctuations af and a^. These can be combined into 63 as discussed above. For simplicity we 
set all other fluctuations to zero and only study 63 here. There is one pure gauge solution 
a = {—iu!,ik), and one non-trivial solution determined by a single parameter, which we 



-19- 



choose to be a^{z = zh) at the horizon z// = 1. Similar to the previous case we define the 
Green's function 

iR _ (vev) _ 63 - \og{z)zd:,e:i 



G 



C3C3 



(source) 



zdze^ 



(4.20) 



Note that satisfy coupled equations and would require a more careful treatment. For 
63 however we numerically find a linear dispersion relation 



UJ 



±k. 



(4.21) 



as seen in figure 4; in formal agreement with [11]. However, note that the authors of [11] 
choose to interpret their correlation functions as current-current correlation functions. In 
contrast to that we interpret our two-point functions according to the alternative setup as 
gauge field correlators. In this context it is reassuring that the gapless mode propagates 
with the speed of light c = 1 at finite temperature and charge density. It behaves just like 
a free photon. 




Figure 4: Linear dispersion relation from our numerics: The quasinormal mode with lowest 
energy of the bulk field E'^ in the normal phase. This corresponds to the lowest-lying pole of the 
corresponding Green's function of the dynamical electric field e = oja^ + kaf. 



Broken phase at vanishing spatial momentum k = Here we are going to demon- 
strate that there is a non-trivial gapless mode appearing in the spectrum of our (1 + 1)- 
dimensional theory. We work at finite W, $ and w, but vanishing momentum k = 0. Then 
we have three coupled fields a}, af, a^, but these can be reduced to one gauge- invariant 
field being dual to d^, and obeying a single boundary condition (apart from the infalling 
condition), as discussed in section 4.2 above. 

We consider the Green's function of the dynamical gauge field combination 0^ in terms 
of the gauge-invariant electric field = coa^ [18] 



G: 



R 



vev 



OJ 



(source) 



OJ 



log{z)zdzal 
zdza% 



(4.22) 
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Due to the presence of a Goldstone mode from the spontaneous breaking of the ^7(1)3, 
we expect that at T = Tc there has to be a gapless mode at a;(A; = 0) = 0. We fohow this 
mode to lower temperatures T < Tc keeping the momentum fixed at A; = 0. The result 
is shown in figure 5. Indeed the mode appears at a; = and moves down the imaginary 
frequency axis as the temperature is lowered. So the Goldstone mode becomes a purely 
dissipative mode at lower temperatures, resembling a diffusive mode.^'^ However, it then 
turns around and approaches a; = again for T — t- 0. Figure 6 shows the temperature- 
dependence of the imaginary gap in this mode. Note that this pole seems to either vanish 
or to become very steep at low temperatures. In the T = OAATc picture of figure 5 it is 
barely visible. 

It may seem surprising that this mode does not acquire a mass gap, i.e. a gap in 
real values of u. This is what one naively expects from spontaneous symmetry-breaking 
in presence of a dynamical gauge field: the Goldstone mode should be eaten by the gauge 
field, making the latter massive. But recall that gauge theories in two spacetime dimensions 
are peculiar as described in section 2.5. Field theory tells us that any non-zero mass for 
the gauge field would also contribute to the chiral anomaly through the relation given 
diagrammatically in figure 1. However, holographically the chiral anomaly is dual to a 
non- vanishing Chern-Simons term of appropriate dimensionality in the bulk action [24-26]. 
Assuming that there is no other way in which the bulk physics can generate an anomaly, 
we conclude that our bulk action (1.1) can not give rise to massive gauge fields. This 
assumption seems safe since the chiral anomaly is a non-perturbative topological quantity 
in field theory. It necessarily needs to be dual to a topological quantity in the dual geometry. 
A pure Yang-Mills term as in (1.1) can not give rise to topologically distinct solutions. 

Figures 7 and 8 show the two-point function of the dynamical field strength e^, or in 
other words the electric field on the boundary. We find an energy gap of size cogap in the 
spectral function, which is presumably related to the energy needed for breaking a charge 
carrier out of the condensate. 

In fact we could employ the picture of a photon propagating through a superconductor 
here. A zero value of that propagator in the frequency-gap means that the photons do 
not propagate. Thus we indeed sec that the J7(l)-interactions become short-ranged here 
although the photons do not become massive. 

Our currents are not dynamical. However, as suggested in [18] we may formally define 
a conductivity via Ohm's law as 

a = i^. (4.23) 

Figure 10 shows this quantity. There is a prominent resonance-like spike near a; = 10 in 
our conventions. From the numerics it is not clear if this spike has finite height. Note the 
surprising match between the quantity uj~^GP- and the Drude model shown in figure 9, 
and discussed in appendix B. We have also tried fitting the formal conductivity given in 
figure 10 using the Drude model but did not succeed. 

Note that our present setup and field configuration are formally similar to the case of 
the holographic s-wave superconductor studied in [43]. But note also the distinct inter- 
pretation: in [43] the fluctuations of gauge fields have been interpreted using the standard 

similar mode was observed in the holographic s-wave superconductor in three spacetime dimen- 
sions [40]. 
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Figure 5: Dynamics in the broken phase: Poles of the correlator Gfs-s (derived from the bulk 
fluctuation a'^) in the complex frequency plane are visible in these contour plots of the spectral 
function oc ImG^s-a. A gapless mode with non-trivial dynamics appears in the broken phase. The 
temperature is lowered from T to T = 0.91Tc, 0.73Tc, 0.57Ta, OAAT^. 
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Figure 6: The imaginary gap: imaginary part of the frequency of the mode which was gapless at 
T^Tc. We keep k = fixed. 



ordinary quantization. In particular this interpretation yielded a superconducting charge 
density which was negative and divergent near T(.. We find exactly the same divergent 
behavior of the corresponding superconducting density from our fluctuations if we use that 
standard interpretation, see figure 12. If we choose the alternative quantization interpreta- 
tion, the formally-defined superconducting density derived from (4.23) via Ug = lim lm{uja) 

uj—>-0 

seems well-behaved at all frequencies, except for one divergence at T* ^ O.bTc, see figure 11. 
This divergence occurs where the graph in figure 12 crossed the real frequency axis. Note 
that the result displayed in figure 11 is essentially the inverse of figure 12. 

This somewhat odd behavior of the formally-defined conductivity suggests that either 
we should not interpret cr as a conductivity, or our (1 -|- l)-dimensional system shows very 
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Figure 7: Spectral function in broken phase shows frequency-gap Wgap ~ 10: Imaginary part of 
the Green's function of the dynamical electric field on the boundary as a function of frequency uj 
at T = 0.169Tc. In the gap the gauge field (our analog of a photon) essentially does not propagate. 
This can be interpreted as the gauge-interaction (our analog of electromagnetism) becoming short- 
ranged. 
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Figure 8: Broken phase: Real part of the Green's function of the dynamical electric field at 

T = o.iegTc. 

peculiar behavior, such as a divergent superconducting density at a particular temperature 
T* = 0.5Tc. 

Broken phase at non-vanishing spatial momentum k It is straightforward to 
find numerical solutions to the full coupled system of equations given in (4.1). However, the 
interpretation of these solutions is not obvious. All the modes are coupled, which means 
that in the dual field theory all the operators mix along the renormalization group flow. 
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Figure 9: Broken phase: Fit of our field strength correlator divided by oj^ to the Drude model. 
The relevant parameters are ctq « 5.8288 x 10^, r w 3.2026 x fO^ when T/T^ w 0.f69. 



Re[o-] 
0.4 r 




5 10 15 20 



Figure 10: Broken phase: The formally-defined conductivity versus frequency at T = O.lGOTc- 

In principle we would like to find quasinormal modes and their dispersion relations 
'^n,QNM(^) in this system. However, we can not simply use the prescription given in [41]. 
This prescription needs to be adjusted to the fact that the sources appear in the logarithmic 
terms of the bulk field expansions near the AdS-boundary. We leave this point for future 
work. 
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Figure 11: The formally-defined superconducting density diverges at T* w 0.5Tc when the 
alternative interpretation is employed. 

5. Discussion 

In this work we have shown that a 
black hole in Yang- Mills theory in 
AdSs can grow vector hair, as seen 
in figure 2. This solution with a 
non-trivial profile for the charged 
gauge field is thermodynamically 
preferred at low temperatures, as 
seen from the lower free energy in 
figure 3. This solution corresponds 
to a boundary theory which de- 
velops a non-zero vacuum expec- 
tation value for a dynamical gauge 
field (21a;) / 0. Therefore the local 
gauge symmetry of this boundary theory is broken spontaneously, leading us to interpret 
this phase as a p-wave superconductor in (1 + 1) dimensions. Formally the symmetry 
breaking in presence of our dynamical gauge field occurs just like in the setups with non- 
dynamical gauge field which yield superfluids. However, we are not able to directly compare 
the two processes within our gravity model since we are restricted to the alternative quan- 
tization, i.e. to the dynamical gauge field case. See section 3.3 for an extended discussion 
of this point. 

Further we have studied the dynamics of this model and have shown that at least one 
non-trivial gapless mode exists in the broken phase. This may allow a low-energy effective 
analysis of this setup in what we defined in section 2.2 as "non-equilibrium field dynamics" 
(this is analogous to the hydrodynamic description in the standard quantization). The 
origin of this gapless mode is the Goldstone boson which appears at the transition tem- 
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Figure 12: The supercounducting density diverges at 
T — Tc when the ordinary interpretation is employed. It 
vanishes at a particular T* « 0.5Tc. 
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peraturc. Naively we would expect this mode to give mass to our U{1)3 gauge field, just 
like the photon becomes massive in a traditional superconductor. However, our numerical 
studies suggest that our gauge field does not acquire a mass, in contrast to our naive ex- 
pectation. Instead, our Goldstone mode below the transition acquires a larger and larger 
damping at lower temperatures, as was seen in figure 6. 

This fact again confirms our interpretation of these two-point functions as correlators 
of dynamical gauge fields (i.e. gauge field propagators) in the following way: We know 
from the Schwinger model that in two spacetime dimensions a C/(l) gauge field can only 
have mass if the theory also has a chiral anomaly, see discussion in section 2.5. However, 
we have forced our anomaly to vanish by forcing the Chern-Simons coupling to vanish in 
the action (1.1). Therefore we conclude that our gauge bosons can not acquire a mass in 
the present setup. See section 4.3 for a detailed discussion. It would be very interesting to 
investigate the action (1.1) after addition of a Clicrn-Simons term. We suspect that this 
will allow a mass generation by a Higgs-like mechanism below the critical temperature.^^ 
In fact, the size of the condensate may be linked to the size of the anomaly, see [36] for 
a similar setup in AdSs. The significance of the formally-defined [18] conductivity (4.23) 
requires more investigation. 

Wc also see the gauge boson propagator (virtually) vanish in the superconducting phase 
at low frequencies, see figure 7. This STiggcsts that the gauge interactions mediated by OTir 
field theory gauge boson di^ (our analog of the photon) become short-ranged. In other 
words, we see that our analog of the electromagnetic interaction becomes short-ranged. In 
this context it would be interesting to perform a calculation at finite momentum, which 
may allow to derive the penetration depth of this superconductor. As a general conclusion 
this holographic superconductor in the alternative formulation shows signatures which are 
strikingly similar to the ones that are known in the ordinary formulation. For example the 
gap in our gauge boson propagator (figure 7) resembles the (conductivity) gap observed in 
the current-current two-point functions in the ordinary formulation, see e.g. [3]. However, 
we stress again, that the interpretation of these signatures and quantities is different in the 
ordinary and the alternative formulations. For example, the gap in the ordinary current- 
current correlators is interpreted as a gap in the (single electron) conductivity. In our 
alternative setup the gap in the corresponding gauge-field two-point function is interpreted 
differently as explained at the beginning of this paragraph. 

Although our action (1.1) preserves parity, our ground state solution spontaneously 
breaks this symmetry. With this property our system may serve as a toy model for un- 
conventional superconductors, which violate parity in their ground state through the order 
parameter. An experimentally well-studied example for such a system in 2 + 1 dimensions 
is the strontium ruthenate compound Sr2RuOi [44]. In this context it may also be inter- 
esting to consider an analog to our setup in AdS4. Note that in the ordinary quantization 
the obvious (p-l-ip)-wave solution which breaks parity in AdS4 is thermodynamically not 
favored according to the analysis in [36]. 

^^In the conclusions of [18] it was noted that for such a theory a perturbative instabihty is ruled out 
due to flatness of the gauge connection. However, non-perturbatively there may exist a thermodynamically 
preferred phase with a (vector) condensate. 
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None of our fluctuations indicates an instability of the p-wave phase. However, a study 
of the full mode spectrum of the coupled system at finite spatial momentum turns out to 
be involved. In particular this requires an extension of the (determinant) method proposed 
in [41] in order to account for the logarithmic terms in the AdSa-asymptotics. This is 
postponed to future work. 

Our paper leaves some other loose ends which are worthwhile being picked up. In 
particular one could repeat our analysis in AdS4, where both quantizations are possible, and 
thus can be compared directly. This should allow to make the relation of Green's functions 
in the alternative quantization to the ones in the ordinary quantization explicit. In AdS4 we 
also have a magnetic field, which can destroy the ordered state, and a dynamical response 
should be observable. So the Meissner-effect could be observed directly. For this purpose 
the authors of [5-8] provide stationary solutions (in the s-wave superconductor), while we 
suggest to examine the dynamics of their backgrounds via the bulk field fluctuations along 
the lines of our present paper. 

In analogy to the standard hydrodynamic study of higher-dimensional setups, sec for 
example [35,45], one could study the alternatively quantized setups at low frequency and 
momentum. Such a study should yield results resembling those of traditional thermal 
quantum field theory, and field dynamics out of equilibrium. Due to the simplicity of our 
setup one may hope for analytic results in the spirit of [46]. 

In fact our investigations suggest that the alternative interpretation may be interesting 
for a wide class of holographic fluids: It should be possible to study each of the holographic 
fluids which have been investigated in AdS4, but now with the Neumann boundary con- 
dition, i.e. in the alternative quantization. Our approach should be applicable to all 
holographic fluids, such as holographic quark-gluon plasma and all superfiuids, and super- 
conductors (s-wave, p-wave, d-wave), which have been found in AdS4. 
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A. Gauge-Invariant and Gauge- Covariant Fields 



A.l Broken Phase 

In the broken phase (superconducting phase), the gauge transformation of the perturba- 
tions are given by 



-iojt+ikx a \ 



(A.l) 



where, a" can not depend on z because we have chosen an axial gauge. Therefore a" has 
constant components. 

One way of writing these gauge transformations of from (A.l) expUcitly is 



( aW ( ik 
a 
a 



a- 



\ 

ik -W{z) 

W{z) ik 

-iuj —$(2;) 

$(z) -iuj 

V -ioj J 



a 



X a" 



a 



(A.2) 



Actually, there are various gauge invariant linear combinations for one typical perturbation. 
For example, we can calculate the linear gauge invariant combinations of in general form 
as follows: 



= o?^ + hal.+ j a1. + ma\ +na1 +pal. 



(A.3) 



where {h, j,m,n,p) are coefficients without a^. Because of Sa^ = 0, we can have 7 types 
of gauge invariant combinations: 

2. h = m = n = 0, p = - ^"-fc ,j =- 



k ' 



3. h = p = 0, m 



(W(^)2-fc2)$(2) 

VK(z)(a;^-*(z)^) 



ik . 
W(z)' 



A.j = m = 0, h = ^^,p=ln = -i 



iW{z) . 



5. J = n = 0, /i 



few 

k^{z) ~ w'"^ 



6. m = p = 0, h = — 



(fc2_H/(2)2)${z) _ i(fc2-W^(2)2) 



n 



ik . 

W{z)' 



^ , tj(fe2_ty(2)2) 

7. n = p = 0, h = ~ IwizMz) ^ "^ 



ujW{z) 

k^-Wjz)^ ■ _ ik 
W{z)<l>{z) Wlzj- 



In fact, one can also construct the linear gauge-invariant combinations for other perturba- 
tions in the same spirit. We will not show them explicitly here. 
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A. 2 Normal Phase 

In the normal phase, only part of the SU{2) is broken and a U{1)^ remains intact. Then 
the fluctuations transform as 

5al = d,X^ + e^'%AlX<^ + a^A^Si) , (A.4) 

where is the gauge transformation parameter of the J7(l)3. From this we obtain 

Sa] = -iuX^ - ^{z)X^ + alls? , 

Sal = -i'^>^^ + ^{z)X^ - a* , 

5al = -iuX^ , 

Sal = ik^^ + a^A^ 1 

Sal = if^^'^ - «x-A^ > 

Sal = ikX^ ■ (A-5) 



Define ef = al± ia^, = =b ial which are in the fundamental representation of U (1)3, 
and is in the adjoint. Using the above transformation, we obtain 



Set = dtX-^ + i^{x)X+ - iA^ef , 
Se^ = dtX" - i^{x)X~ + lA^e^ , 
Se+ = d^X+ - iA^e+ , 
Se~ = dxX~ + iA^e~ , 

Sal = • (^-6) 

Prom the af and a^, we can construct a i7(l)3-invariant combination with the help of 
Fourier transformation (e~*'*'*+*'^^/(t, x) f{oJ,k) as before), that is 

63 = kal + ujal , (A.7) 

which leads to 5e-^ = —ikcoX^ + ikujX'^ = 0. 

On the other hand, we can build C/(l)3 gauge-covariant combinations for ef and e^, 

e+ = {^{z)-u)et-ket, (A.8) 
el = {^iz)+cj)e- +ket. (A.9) 



This leads to 



and 



Sel = ($(z) - Lo)5e+ - kSej = -iA^el , (A.IO) 



Se^ = {^{z) + uj)Se~ + kSe^ = iA^el . (A.ll) 

With these gauge-covariant fields and and the gauge-invariant field 63 in the normal 
phase, we can rewrite the on-shell action. From the on-shcll action we can see that 63 
decouples from all other fields while and ej^ are coupled to each other. 
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B. Ordinary Interpretation 



B.l Drude Model 

The quantity is identical to the conductivity one would define in the ordinary 

interpretation of our setup. Therefore, we tried fitting this using the Drude model 

Re(c7Drude) = -, , ^"o 2 " (^-1) 
1 + UJ'^T'^ 

where, gq = nq^r/m is the DC conductivity of axx, while n, q, m, r are respectively the 
electron's number density, charge, mass and the mean free time between collisions. As seen 
from figure 9 this works surprisingly well at low frequencies. We have no good explanation 
for this agreement. 

B.2 Superconducting Density 

Here we assume that the ordinary interpretation of sources as fixed boundary gauge fields 
is correct. In this case our two-point functions would be correlators of dynamical currents. 
In particular this would allow us to compute the superconducting density using the formula 

Us = lim Re G.f (w; fc = 0) . (B.2) 

Here the Green's functions is obtained from the infalling solution of the bulk field given 
in (4.9). The result is shown in figure 12. We observe the same negative divergent behavior 
as seen in [43]. We take this as a numerical evidence for the fact that the alternative 
interpretation should be applied to this kind of setup. See section 4 for a detailed discussion 
of the alternative interpretation. 
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